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The critical behaviour of a Random Fiber Bundle Model with mixed uniform distribution of 
threshold strengths and global load sharing rule is studied with a special emphasis on the nature 
of distribution of avalanches for different parameters of the distribution. The discontinuity in the 
threshold strength distribution of fibers non-trivially modifies the critical stress as well as puts a 
restriction on the allowed values of parameters for which the recursive dynamics approach holds 
good. The discontinuity leads to a non-universal behaviour in the avalanche size distribution for 
smaller values of avalanche size. We observe that apart from the mean field behaviour for larger 
avalanches, a new behaviour for smaller avalanche size is observed as a critical threshold distribution 
is approached. The phenomenological understanding of the above result is provided using the exact 
analytical result for the avalanche size distribution. Most interestingly, the prominence of non- 
universal behaviour in avalanche size distribution depends on the system parameters. 

PACS numbers: 64.60.Ht, 81.05.Ni, 46.50.+a, 62.20.Mk 



I. INTRODUCTION 

Breakdown phenomena in nature has captured the at- 
tention of scientists for years A study of this phenom- 
ena plays a major role for the prediction of failure and 
design of materials and structures. One of the paradig- 
matic model mimicking the fracture processes is random 
fiber bundle model (RFBM) which is simple yet subtle 
enough to capture the essential physics of the breakdown 
phenomena. 

Random fiber bundle models @, d, 0, @ have been 
studied extensively in recent years. Typically a RFBM 
consists of N parallel fibers with randomly distributed 
threshold strength {ath) taken from a given distribution. 
If the stress generated due to an external force is greater 
than <j t h of a fiber, it breaks. The dynamics of the model 
is initiated by applying a small external force just enough 
to break the weakest fiber present in the bundle. The 
load carried by this broken fiber is shared amongst the 
remaining intact fibers following a load sharing rule caus- 
ing further failures. When no further failure takes place, 
the external force is once again increased quasistatically 
to break the weakest intact element present in the bun- 
dle and the process continues till the bundle breaks down 
completely at an external stress called the critical stress. 
Even though the threshold distribution of real materials 
may not be known exactly, in theoretical models the dis- 
tributions are usually approximated by either a uniform 
distribution or a Weibull distribution Q . 

The avalanche size is defined as the number of bro- 
ken fibers between two successive loadings. The distri- 
bution of avalanche size turns out to be a key factor in 
charaterising any breakdown phenomena. Hemmer and 



Hansen[6| studied the avalanche size distribution D(A) 
of an avalanche of size A in a RFBM with the global load 
sharing (GLS) scheme, in which the additional stress due 
to a broken fiber is distributed equally to the remaining 
intact fibers. They established a universal power-law dis- 
tribution in the large A limit given as D(A) ex A - ^ with 
£ = 5/2. 
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FIG. 1: Mixed Uniform Distribution 

RFBM with threshold strengths which are continuous 
and are uniformly distributed between to o\ and also 
between ai to 1 [3|, |4[ , have been studied separately. But 
what happens when two such bundles are merged is not 
known, especially when ax < a% such that there exists a 
discontinuity in the threshold strength distribution. In 
this paper, we investigate the role of such a discontinuity 
on the critical behaviour of a RFBM. The distribution of 
threshold strength of fibers used in the present work is 
given as (See Fig. 1) 
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The discontinuity as denned above introduces dilution 
in the model in the sense that two types of fibers sepa- 
rated by a gap in their threshold distributions coexist in 
the same bundle. It is shown below that this disconti- 
nuity plays a crucial role in the dynamics of the model. 
Here, a fraction / of fibers belong to the weaker thresh- 
old distribution with strengths uniformly lying between 

and CTi (Class A) whereas the remaining fraction of 
fibers have stronger threshold strengths between (72 and 

1 (Class B). Clearly, (cr 2 — ax) is the measure of the dis- 
continuity which vanishes in the limit of purely uniform 
distribution [J]. 

In a recent paper, Pradhan, Hansen and Hemmer [7J 
showed that for a bundle which is close to the complete 
break down (i.e., imminent failure), a crossover in £ from 
a value 5/2 to 3/2 is observed when the threshold distri- 
bution approaches the critical distribution. Critical dis- 
tribution in their case is the distribution in which the 
lowest threshold of the remaining intact fibers is equal 
to half that of the strongest. This crossover has also 
been observed with other load sharing rules Q. We are 
however interested in looking at the distribution of to- 
tal avalanche size and a crossover from a power-law be- 
haviour with £ = 5/2 to a non-mean field, non- universal 
behaviour for smaller A is observed near a critical distri- 
bution. This interesting behaviour is due to the presence 
of class A fibers. In the present model, the critical dis- 
tribution corresponds to 02 = 0.5. We emphasize that in 
Ref. [7] , the proximity of xq (the threshold of the weakest 
intact fiber) to the the critical threshold (=0.5) leads to 
the crossover behaviour whereas in the present case it is 
the proximity of 02 to the critical value (0.5) which is at 
the root of the observed crossover behaviour. In a sense, 
(T2 is playing a role analogous to xq in our model. 

The paper is organized as follows: We have already 
introduced the model above. The critical stress and ex- 
ponents are obtained in Sect. 11(A) with a special em- 
phasis on the avalanche size distribution in Sect. H(B). 
The concluding remarks are presented in Sect. III. 



II. RESULTS AND DISCUSSIONS 

A. Critical stress and exponents 

To study the dynamics of failure of fibers, we use the 
recursive dynamics approach Q. If a fraction / of the 
total fibers belong to Class A and the remaining 1 — / to 
class B, then the uniformity of the distribution demands 

1 f ai 
f = T— 7 s / da; 
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The above equation provides a relationship between /, 
<7i and (72 and at the same time puts a restriction on the 



allowed values of the parameter ax as shown below. Any 
value of (7i > / leads to a value of 02 smaller than o\ 
which is not an acceptable distribution (see Eq. 1). 

We now define Ut as the fraction of unbroken fibers 
after a time step t. Then the redistributed stress at the 
instant t is a t — F/N t — a/U t where the applied force 
F = Na and N t — NUt- The recurrence relations be- 
tween U t , Ut+i and between a t , er t+1 for the GLS are ob- 
tained as 



U+i = 
and <Jt+i = 



1 - P{°t) 
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where P{a t ) is the fraction of broken fibers with the re- 
distributed stress at, and is given as 



p{ath)d(Jth- 



This dynamics propagates until no further breaking takes 
place. It should be emphasised that the initial load is so 
small that the redistributed stress is always less than (72 
and therefore fibers from class B cannot fail. Thus to 
initiate the breaking of class B fibers, the redistributed 
stress at a later time t must exceed 02. The fixed point 
solution for U (= U*) and a(— a*) at which no further 
failure takes place can be obtained using the standard 
technique of solving the above recursive relations. By 
substituting P(a t ) in Eq. (3) we get 



U 



t+i 



P(~) 



1 



l-(<7 2 -<7l) l-(<7 2 -<7l) U t 

At the fixed point, 
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giving the following stable fixed point solutions 



U* 



1 



2(1 - (a 2 - ax)) 1 



1 ] and 

C 



1 1 / a 

2 ~ 2V 1 ~ ~a c 

1 

4[l-((7 2 -(7l)]' 



with 



(4) 
(5) 



If the external applied stress is less than a c , the system 
reaches a fixed point. For a > a c , the bundle breaks 
down completely as both the U* and a* become imagi- 
nary. Equation (4) suggests that the redistributed stress 
a* attains the maximum value (=0.5) at a c . The critical 
stress of the mixed model varies with the gap [pi — ax) 
and reduces to the value a c = 1/4 for the uniform dis- 
tribution as (02 — 0i) — » 0. In deriving Eq. (4), a t is as- 
sumed to be greater than (72 while the applied stress when 
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plotted against the redistributed stress shows a disconti- 
nuity at <j\. Beyond <j\, the external force is increased 
to break the fiber with threshold strength 02, i-e., the 
gap in the threshold distribution is also reflected in the 
constitutive behaviour of the model. 

The discontinuity on the other hand imposes some re- 
strictions on the parameters for this calculation to be 
valid. Since the maximum value of the redistributed 
stress is equal to 0.5, cr 2 must be less than 0.5 so that 
some fibers from class B also fail at the critical point. 
The condition a% < 0.5, eventually restricts the value 
of critical stress to be less than 0.5 for any chosen dis- 
tribution. However, 02 — 0.5 is a limiting case when 
the redistributed stress at the critical point marginally 
reaches class B fibers. In short, we must have 

<J\ < f,<J2 < 0.5 . 

One can define an order parameter O associated with the 
transition [4| as shown below: 



O = 2[1 - (a 2 - ax)]U* - 1 = (a c - a) 1 ' 2 = (a c - af. 

The order parameter goes to as a — * a c following a 
power law (a c — er) 1 / 2 . Susceptibility can be defined 
as the increment in the number of broken fibers for an 
infinitesimal increase of load. Therefore, 



dm 

X = —r- where 
da 



= N[l - U*(a) 



Hence, 

X oc (cr c - er)~ 1/2 = (a c - er)~ 7 . 

The exponents f3 and 7 stick to their mean field values 
0, [1[ i.e. (3 = 7 = 1/2. We conclude that though the 
discontinuity alters the critical stress, the critical expo- 
nents remain unaltered. It is to be noted that the model 
reduces to the already obtained results in various lim- 
its. For example, if class A fibers are absent (/=0), 
a c = 1/4(1 — 02) an d an elastic to plastic deformation 
is observed [4|. 

B. Avalanche size distribution 

Let us now focus on the avalanche size distribution 
exponent £. Below is shown some of the allowed distri- 
butions which satisfy the above mentioned restrictions. 
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0.30 


0.25 
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0.30 


0.29 
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0.26 
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0.40 


0.35 


0.47 


0.29 



We study the avalanche size distribution numerically 
by using the method of breaking of the weakest fiber [g]. 
In the simulation, the fibers are arranged in an increasing 
order of their threshold strengths. An external force suf- 
ficient to break the weakest fiber is applied and the load 
due to the breaking of this fiber gets redistributed among 
the remaining intact fibers following the GLS scheme. 
The number of failed fibers for a fixed external load is 
recorded till the dynamics reaches a fixed point. There- 
after, the external load is increased further and the above 
process is repeated till the critical stress is reached. 

Following interesting observations (Fig. 2) are clearly 
highlighted, (i) For the cases 1, 2 and 5 of the table, 
the avalanche size exponent is 5/2. (ii) For the cases 
3,4 and 6, there is an apparent power law behaviour for 
smaller A with the exponent which is found to depend 
on the system parameters. In the examples chosen here 
the exponent happens to be close to 3. For larger A, 
however, we retrieve the universal mean field behaviour 
with £ = 5/2. Also, (hi) an increase in the region with 
£ w 3 is seen as 02 — * 0.5. These observations establish 
the following: (i) there is a non-universal behaviour of 
D(A) in the small A limit, (ii) there is a crossover to 
the universal behaviour in large A limit, and (iii) the 
crossover behaviour is prominent as 02 — > 0.5. 

The above mentioned results can be explained by ex- 
tending the analytical result for the avalanche size distri- 
bution obtained by Hemmer and Hansen Q to the mixed 
model. The general expression for the avalanche size dis- 
tribution with GLS is given as 
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where x is the redistributed stress and Q(x) is the frac- 
tion of unbroken fibers at x. The upper limit of the 
integration (x c ) is the redistributed stress at the criti- 
cal point. The right hand side of Eq. (6) is broken into 
two parts, -Di(A) and D 2 (A) for the mixed model with 
any allowed values of o\ and o<i (where p(x) = p = 
1/(1 - cr 2 + o-i)): 

As long as the redistributed stress is restricted to the 
class A fibers, 



Q{x) = 
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When the redistributed stress belongs to the second block 
(class B), we have 
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FIG. 2: Fig. (2a) corresponds to easel of the table where 
no crossover is observed. (2b) shows a crossover from a non- 
universal behaviour to the universal behaviour with £ = 5/2 
as A is increased (case4). Fig. (2c) shows an increase in region 
with non-universal behaviour (case6). A line with slope -5/2 
is drawn in each figure for comparison. In Fig. (2c), a line 
with slope -3 is also shown to indicate the apparent power-law 
behaviour. Fig. (2d) compares the simulation and numerical 
integration data for case4 which match identically with each 
other. 
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A close inspection shows that Eq. (8) corresponds to a 
situation where class A fibers are absent and the results 
match exactly with Pradhan et.al 0] where a crossover 
in the avalanche size exponent from 5/2 to a value 3/2 is 
observed as o-i — > 0.5. Here, the presence of class A fibers 
necessitates the simultaneous study of Eqs. (7) and (8) 
for the total avalanche size distribution. For distributions 
which are closer to the critical distribution, Eq. (8) yields 
exponent £ = 3/2 for smaller A, but the contribution 
from the Eq. (7) modifies this small A behaviour. 

The results obtained by numerically integrating 
Eqs. (7) and (8) match identically with the numerical 
simulation results [Fig. (2d)]. We observe a sharp fall of 
D\ (A) [Eq. 7] which points to the fact that the avalanches 
of smaller sizes are contributed mainly by the breaking 
of class A fibers (Fig. 3). The point of intersection of the 
two integrals yields A c which is the value of the avalanche 
size where the crossover takes place. The variation of A c 



FIG. 3: Figs, a, b and c present the numerical integration 
results showing an increase in A c (defined in the text) as &i — > 
0.5. The symbol corresponds to log 10 D( A) /N=log 10 (A (A) + 
D2(A)) whereas log 10 Z?i(A) and log 10 D2(A) are drawn to 
compare their relative magnitudes. 



for three different cases are shown in Fig. 3. Incidentally 
in the examples presented here, A c increases with o^- 
If we look at a particular critical distribution (Fig. 4), a 
crossover from an apparent power-law behaviour of D(A) 
with exponent £ close to 3 for small A to a universal be- 
haviour with £ =3/2 is observed. 

We now discuss the numerical integration and simula- 
tion results in the light of Eqs. (7) and (8). With a change 
of variable x —* x/(l — <T2 + <xi — x) , we can rewrite Eq. (7) 
in the following form 
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where x m — <xi/(l — 02) ■ The argument of the expo- 
nential term has a maximum at x = 1 which is outside 
the range of integration and hence the saddle point in- 
tegration method can not be applied [9J ■ Expanding the 
term l/(l + a;) 2 ina power series and then using the in- 
complete gamma function (lp|. we arrive at the following 
result (see the Appendix) 
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FIG. 4: Avalanche size distribution for a particular critical 
threshold distribution with / = 0.4. Thick line corresponds 
to a slope of -3 and dotted line to -3/2. The asymptotic 
behaviour clearly is D(A) ~ A _3//2 . 
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(10) 



The leading behaviour of the infinite series (10) is 



A- 5 /2 e (i- 



which justifies the non-universality 



observed in numerical simulations in the small A limit. 
The question therefore remains why does the non- 
universal behaviour become prominent as a% — > 0.5. The 
behaviour of D 2 (A) [Eq. (8)] for smaller A(<< A c ) is 
at the root of this. When a 2 — > 0.5, Z?2(A) goes as 
A -3 / 2 and hence the contribution from Di(A) wins over 
to produce a prominent non-universal behaviour. The 
small A behaviour of D(A) is therefore non-mean field 
and non-universal when 02 — * 0.5. On the other hand if 
02 << 0.5, x m is relatively smaller and the contribution 
from class A fibers decays very fast as A increases and 
one observes a mean-field universal behaviour almost for 
the entire range of A. We also observe a non-powerlaw 
behaviour for small A when o\ is very small and &2 is 
close to criticality as is expected from the analytical re- 
sult. However for a given cr 2 (~ 0.5), larger x m leads to 
a non-universal behaviour up to larger A. The crossover 
value A c is thus roughly given by the value of A for which 
£>2(A) crosses over from a A _3//2 behaviour to A -5 / 2 be- 
haviour. 



III. CONCLUSIONS 

In conclusion, we have studied a mixed fiber bun- 
dle with a discontinuous but uniform threshold distri- 
bution and GLS. Discontinuity leads to a functional de- 
pendence of the critical stress on the system parame- 
ters a\ , er 2 and / and also imposes restrictions on the 
allowed values of these parameters. Although the critical 



exponents are unchanged, there is a non-trivial change 
in the burst avalanche distribution behaviour where dis- 
continuity leads to a non-universal, non-mean field be- 
haviour for small A. We would like to emphasise that 
non-universality becomes prominent only when 01 — » 0.5. 
For large A limit, the behaviour is however universal and 
mean field. If / = or 02 << 0.5, the non-universal be- 
haviour completely disappears. The non-universality in 
D(A) is also seen for other distributions @. The beauty 
of our model is that the non-universal behaviour is tun- 
able with the system parameters and there is a crossover 
from non-universality to universality in the limit of large 
A. One should also note that the imminent failure of 
the bundle (i.e., final stages of the breakdown process) is 
the same as in Ref. Q because the effect of class A fibers 
essentially vanishes in that limit. 
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Appendix 

In this appendix, we shall indicate how to arrive at 
Eq. (10) of the text starting from Eq. (9). Equation (9) 
can be written as 
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The integral can be further written as 
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Using the power series expansion of Incomplete Gamma 
Function 
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Rearranging terms, we get x m (x m A) 



x , Q (A + q + 1)(A + q + 2) . . . (A + q + 1 + k) ' 

/(A) = e-*™ A x££(-l)*(<?+l)a4x 

9=0 

(Yt-, 77^ / a TT i.e., Eq. (10). 
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